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Abstract: We study thermal corrections to the dispersion relations of massive funda- 
mental particles immersed in weakly coupled non-Abelian plasmas. The cases covered 
include quarks in the QCD (quark-gluon) plasma, as well as A/" = 2 quarks and scalars 
in an A/" = 4 Super Yang-Mills plasma. We perform the calculations to leading order 
in a weak coupling expansion, and consider all mass scales of the fundamental fields, 
ranging from massless particles all the way to bare masses parametrically larger than 
the temperature. 
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1. Introduction 



The goal of the paper at hand is to study weakly coupled deconfined gauge theories 
at finite temperature, and in particular to determine the leading order thermal correc- 
tions to the dispersion relations of charged fundamental representation particles. The 
dispersion relation provides the on-shell relationship between an excitation's energy 
and momentum, as well as its lifetime, and hence encodes valuable information about 
the nature of the quanta making up the theory. In this paper, we consider the weakly 
coupled limits of both QCD (with Ni flavors or massless quarks) as well as A/" = 4 Super 
Yang-Mills theory (SYM). To these theories, we add one fundamental representation 
test particle with an arbitrary mass, which for QCD is a quark and for SYM an A/" = 2 
hypermultiplet, containing one Dirac fermion and two complex scalars of the same 
mass. Our presentation will cover all mass scales of the fundamental particle, from the 
massless limit to flelds with masses parametrically larger than the temperature. 

The original motivation for our work arose from various studies of strongly cou- 
pled J\f = 4 SYM, where the dispersion relations of fundamental particles had been 
determined using the AdS/CFT correspondence (see e.g. Ref. In the course of our 
project, we however noted the surprising absence of some of the corresponding results 
in the QCD literature — in particular in the limit of heavy quarks, with M > T — and 
therefore decided to include this theory in our presentation as well. Our paper thus 
contains a self-consistent study of the leading order (LO) thermal corrections to the 
dispersion relations of fundamental particles in both theories and for all mass scales. 
We note, however, that a subset of the results presented in Section 5 on QCD can be 
found from the literature, in particular in Refs. P, |], ^. In addition, in the massless 
limit of QCD, even the NLO thermal corrections to the dispersion relation of a static 
quark have been recently determined 0. 

For QCD with massless and light quarks at soft {0{gT)) momenta, our flndings 
confirm the earlier results of Refs. P, ^, which we furthermore observe to be closely 
analogous to the ones we obtain for fundamental quarks in A/" = 4 SYM. For the 
latter, one merely needs to change the value of the soft mass parameter appearing in 
the QCD result, but its form stays otherwise intact. Moving then on to SYM with 
either bare masses or three-momenta of order T, we find a simple dispersion relation 
E = \/ 'p^ ^ + (JM^ for both quarks and scalars, of which the latter in fact obey this 
relation even in the soft regime. For SYM, the thermal mass shift (5M^ is independent of 
the momentum p, whereas in QCD, the dispersion relation is seen to have a non-trivial 
form for practically all bare masses and momenta. 

In the latter part of our paper, we specialize to the case where the bare mass of 
the fundamental particle is parametrically larger than the temperature, M > T jg. 
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This is a particularly interesting regime to study, as it allows us to compare the weak 
coupling results obtained in SYM to the corresponding strong coupling calculations 
of Ref. Regardless of the strength of the coupling, we expect to find finite cor- 
rections to the dispersion relation due to infrared (IR) physics, as thermal screening 
and damping effects of long wavelength degrees of freedom must change the energy 
stored in the fields sourced by the heavy fundamental particle. At weak coupling, the 
IR sensitivity of the dynamics requires us to use resummed perturbation theory in 
determining the dispersion relation, but we argue that the corresponding corrections 
can just as well be computed using macroscopic classical field theory, where thermal 
effects are taken into account through temperature dependent screening masses as well 
as permittivity and permeability tensors. To demonstrate this, we compute the leading 
order non-relativistic dispersion relation both using classical field theory and resummed 
perturbation theory. The two approaches are seen to lead to the same outcome, which 
is later compared to the strong coupling result. 

The paper is organized as follows. In Section 2, we introduce our notation and 
conventions, write down the Lagrangians of the two theories, and explain in detail, 
how we intend to approach the task at hand. Sections 3, 4 and 5 on the other hand 
contain the bulk of our calculations, the two former dealing with A/" = 4 SYM theory 
and the last one with QCD. In these Sections, we assume the mass of the fundamental 
particle to be at most of the order of the temperature, while the non-relativistic limit 
of M ^ T is left to Sec. 6. In Section 7, we then discuss the main implications of our 
results, as well as briefly draw conclusions. While most of the details of our calculations 
are presented in the main body of the text, we have decided to separate the lengthy 
computation of one specific integral, contributing to the self energy of an infinitely 
massive quark, into Appendix A. The treatment in the main text is however entirely 
self-consistent, as we present the entire classical calculation in Section 6. In addition, 
we list the results for some straightforwardly computable sum-integrals in Appendix B, 
but omit the derivations that can be found from several textbooks (see e.g. Ref. 0). 

2. Setup 

In field theory, the dispersion relation of an excitation is obtained by studying the 
location of the pole(s) of the corresponding Minkowski space two-point function. The 
location of the pole encodes information about both the kinematics and lifetime of the 
excitation, thus providing insight into its real time properties. It is, however, both 
possible and convenient to obtain the dispersion relation from time-ordered Euclidean 
correlation functions. As can be easily demonstrated by decomposing them in terms of 
spectral densities, time-ordered Euclidean two-point functions contain the same infor- 
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mation as their real time counterparts, and in particular, upon analytic continuation 
to Minkowski space, have poles in the same locations P, 0. 

The time-ordered Euclidean correlation functions, as well as their poles, in general 
have both real and imaginary parts, but one often finds that the imaginary parts are 
suppressed by higher powers of the coupling constant than the real ones. As the real 
part of the dispersion relation gives the on-shell relationship between the corresponding 
excitation's energy and momentum and the imaginary part its decay width, the fact 
that we only work to the leading perturbative order in the present paper implies that 
we will not be able to extract information on the widths. 

Being solely interested in the finite temperature corrections to dispersion relations 
in this work, we always regularize our sum-integrals by subtracting away their (often 
divergent) zero temperature pieces. Due to only working to leading order in a weak 
coupling expansion, this simply amounts to defining the tree level mass we use in our 
equations as one already containing all zero-temperature quantum corrections, i.e. be- 
ing the renormalized, physical mass parameter. Beyond this, we never need to worry 
about the details of renormalization. 

2.1 Notation and conventions 

Our notation is as follows. We use capital letters P to denote four-vectors, lower case 
bold ones p to denote three-vectors, and the usual lower case letters p to denote the 
magnitude of the latter. Sum-integrals will, as usual, be written as (d = 4 — 2e) 



and their finite parts, where the T = contribution has been subtracted off and d 
subsequently set to 4, as 



Here, the subscript P and the sign + refer to bosonic and {P} as well as — to fermionic 
Matsubara sums, with po = 2mTT and po = {2n + 1)ttT, n G Z, respectively. Unless 
otherwise stated, also the zero components of the external momenta in the graphs we 
consider will be assumed to take these values. The parameter A denotes here the MS 
renormalization scheme scale parameter, but will never show up in the final results. 

The analytic continuation to the Minkowski space is defined by po ^ iS + S, where 
S is an infinitesimal positive real number and S the Minkowski space energy. We denote 
the unperturbed dispersion relation by E{p) = \/p^ ^ M^, and from Section 3 onwards 




(2.1) 
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will work exclusively with the Euclidean metric, gfj_i, = 6fj_u- Finally, a group theory 
constant that will be in frequent use in the following is the Casimir of the fundamental 
representation of SU(Xc)) 

C^5,, = (T»T% = ^^S^r (2.3) 

2.2 The Lagrangians 

The two theories we study in this paper are J\f = 4 SYM theory coupled to one 
fundamental M = 2 hypermultiplet with mass M, as well as QCD with Ni massless 
quarks and one test quark with mass M} They are defined through the following four- 
dimensional Lagrangians, which we give in Minkowski space with — h ++ signature. 

The field content of A/" = 4 SYM consists of one gauge field A^, four Majorana 
fermions ipi and three complex scalars 0p, while the additional M = 2 multiplet is 
composed of two complex scalars $„ and one Dirac fermion uj. All A/" = 4 fields 
transform under the adjoint representation of the gauge group SU(iVc) , while the Af = 2 
sector transforms under the fundamental representation. The A/" = 4 fields are massless, 
and the Af = 4 theory conformal, while the M = 2 fields each have mass M. 

Following Ref. 0, we define 0p = 1/a/2 (Xp + il^), with Xp and Yp hermitian, 
which allows us to write our Lagrangian in the form 

C = Cq + Ci + C2, (2.4) 

with 

£0 = -tr I^F^^F^^ + i,,Ipi,, + (DX,)' + (DF,)' } 

- mi{-D^ + M2)$„ - uj{U) + M)uj, (2.5) 
= tr I - iiJia^^-[Xp, + V^i75/?f^.[Fp, - u) {Yi - 275X1) u 

+ 2V2Im(^ - cuP+^i^i - {^2)HiP+0J + {^i)^2P+uJ - 0jP+i^2^^ 

- 2M($)J,Fi$„, (2.6) 

C2/g' = -\tT {t[XA,XB])' + (-l)"($n)^ ([02,0^] + [03,0^]) $n 

- 4Re(($i)t[02,03]$2) -i|(-l)"($n)^ta<fnr-2|(<l'2)^ta$ir 

- ($„)t{0i,0t}$„. (2.7) 



^Due to the low perturbative order we are working in, the parameter Ni will only be visible through 
its effects on the Debye mass mn of QCD. 
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^^^^^A 

Figure 1: The one- loop diagrams contributing to the self energy of a fundamental a) scalar 
and b) fermion in the = 4 SYM theory, and c) a quark in QCD. The thick and thin dotted 
lines stand for the N = 2 and = 4 scalars, respectively, and the thick and thin solid lines 
for the fundamental and = 4 fermions. The wiggly line on the other hand denotes a gluon. 




Here, D denotes covariant derivatives in the appropriate representations of SU(A^c), 
X = (-^1, Yi, ^2, ^2, ^3, ^3), and a sum over repeated indices is implied. The matrices 
and are given by 



ia2j \-icr2 J V-o-Q 



and they satisfy the algebra 



[a^/?^]=0. (2.9) 
Finally, for QCD the bare Lagrange density reads 

jCqcd = —iiF^.F^'' -^f{lp + Mf)^f, (2.10) 

where F^^ is the field strength, ipf are fundamental Dirac fermions and D is the gauge- 
covariant derivative. The flavor sum over / includes Ni massless flavors with Mj = 
and one test quark with Mf = M. 

3. Fundamental scalars in = 4 SYM 

Let us begin by considering the dispersion relation of a fundamental Af = 2 scalar in 
the Af = A SYM theory. To this end, we note that the Euclidean scalar self energy is 
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given by the sum of the graphs of Fig. |l].a and reads in the Feynman gauge 



2 „ ^ rf^fc / 1 , 1 1 P' + M' 



n(P) - Ag Cf% ^ + ^2 + M2 - 2{P + K)^ + M^l0' ^^"^^ 
, ^ d^k f 1 1 P^ + M^ 1 



{2tiY\K^ K^ + M^ {P + K)^ + M^K\ 
The equation for the pole in the scalar propagator is then simply given by 

+ + n(P) = 0, (3.2) 

from which we wish to solve po- If n(P) is regular in the unperturbed on shell limit 
P^ = — M^, i.e. does not diverge there, then Eq. ( |3.2| ) may be solved to leading order 
using this relation in the evaluation of n(P). With this in mind, we now consider 
the terms proportional to P^ + in Eq. ( p.l|) , which contain the sum-integrals of 
Eq. dRlD of Appendix B. 

To see how these functions behave in the unperturbed on shell limit, we analytically 
continue po to real energies E{p) = a/p^ + by setting po = iE{p) + 5, with 5 > . 
The logarithm appearing in the first term of Eq. (|B.1|) now becomes 

( {E{p)+p)k + t{E{p)-k)5 {E{p)-p)k-t{E{p) + k)5 \ 
°^ V iE{p) + P)k - i{E{p) + k)5 {E{p) - p)k + i{E{p) -k)6j' ^ ' 

which in fact vanishes in the 5 — > limit for all non-zero M, and thus also in the 
M ^ limit. At the same time, the log in the second term of Eq. ( |B.1| ) obtains the 
form 

( M^-E{k)E{p)+pk M^ + E{k)E{p)+pk \ _ ( {k-pf \ 
°^\M^-E{k)E{p)-pk M'^ + E{k)E{p)-pk) ^^\{k + py)' 

which obviously contains an integrable logarithmic singularity at k = p but is otherwise 
regular. It follows that the integrals in Eq. ( |3.1| ) are regular in the unperturbed on shell 
limit, and that the terms in n(P) proportional to P'^ + M^ may be thus neglected when 
evaluating the function. 

Throwing now away the terms discussed above and using Eq. (p.2|) from Appendix 
B to evaluate the remaining sum-integrals, we arrive at the p independent result for 
the self energy, 

n(P^ = -M^) = ^ / dkk-^+ dk4--§^^\ (3.5) 



7r2 



e2/3fc - 1 ' 7o E{k) e^/^^C^) - 1 



6Mt. 
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Figure 2: The form of the leading order thermal correction to the mass squared of a funda- 
mental J\f = 2 scalar field, 6Mg, as defined in Eq. (3^). We give the result as a function of 
M/T, normalized by g^CpT"^. 

This implies that for all values of the bare mass and three- momentum, the scalar 
dispersion relation has the free form 

£ = + M2 + (3.6) 

where only the value of the mass shift SM^ has any dependence on the temperature. 
If the bare mass M is very large in comparison with T, the second integral in Eq. ( p.5|) 
becomes exponentially small and SM^ = ^g'^CpT'^, while in the massless case, the 
two integrals contribute equally to yield a result twice as large. The full function 
SM'^{M/T) interpolating between these two cases is difficult to obtain analytically, but 
is numerically readily available and is plotted in Fig. |[ 

4. Fundamental quarks in A/^ = 4 SYM 

Next, we inspect the case of fundamental M = 2 quarks in the SYM theory, and thus 
consider the leading order fermion self energy as given by the three graphs of Fig. |^.b. 
Continuing to work in the Feynman gauge, we observe that their sum can be written 
as 

d^k i{f + lf[)+M 1 



S(P) = Ag'Cp 



(27r)4 (P + KY + M2 K'^ 



^ ^ d^k ilf[ 1 
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from which we obtain the equation for the pole of the propagator, 

det (z/' + M + S(P)) = 0. (4.2) 
Looking at Eq. ([4.1|) , we note that S(P) may be written in the form 

S(P) = {if + M) a{P) + «f (P), (4.3) 
where the functions a(P) and V^{P) are defined by 



d^k 1 



rf^A; 1 



(4.5) 



(27r)4 (P + 7^)2 + M2 J • 

This imphes that the fermion dispersion relation can be solved from the equation 

(P(l + a(P)) + V{P)f + M^{1 + a{P)f = 0, (4.6) 

which can be further simplified by noting that a simple power counting exercise shows 
that both a(P) and V"^(P) are regular and at most of order g'^T/pQ, when P is on shell. 
The terms in the dispersion equation involving a(P) are seen to be either proportional 
to P^ + M^, and thus vanish at this order, or of the form aP ■ V, which is always 
parametrically small. We conclude that we may thus neglect a(P) altogether in our 
study and write Eq. (|4.6| ) in the form 

{P + V{P)f + = 0, (4.7) 

which motivates us to next analyze the behavior of the functions P • V{P) and V{PY 
in various kinematic regimes. 

4.1 The sum-integrals 
4.1.1 P-V 

To inspect the integral P ■ l^, we first separate the vector V{P) into its bosonic and 
fermionic pieces by writing V = V'^ + V~ , where the latter two functions read 

V-(P) - ««^C.|f (4.8) 
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Completing squares and utilizing the results of Appendix P, we obtain from here 
P-V^ = (4.9) 



27r2 



+ Re 



Re 



dk 



dkk 



m 



2 + 



P^ + M^ 



2pk 



log 



P2 + m2 + 2ikpo - 2pk 
P^ + M^ + 2ikpo + 2pk 

P"^ + + 2iE{k)po - 2pk 



.p2 + M2 + 2iE{k)po + 2pkJ J e^^^^) ± 1 



As discussed already in Section ^ upon taking the real parts of the logarithms, the 
integrals in the terms proportional to P^ + are regular in the limit P^ — M^, 
and may consequently be ignored. This means that in the unperturbed on shell limit, 
we have to leading order in the coupling 



P-V^\ 



P2+M2=0 



dkk 



1 



oo ^2 

dk 







E{k) e/^^(^') ± 1 



(4.10) 



implying that however small the three-momentum p becomes, P ■ are always con- 
stants of order g^. 

4.1.2 V"^ 

In evaluating the V(P)^ term in Eq. ( ^4.7| ), we note that we may from the beginning 
take both p and M to be of order gT or higher, as for larger momenta and bare masses 
this term is certainly subdominant with respect to P- V^. Concentrating first on Vo(P), 
we obtain 



Im 



+ 



2-k'^P 
E{k) + ipo 

m 



dkk 



log 



P2 + M2 + 2ikpo + 2pk 
P^ + M^ + 2ikpo - 2pk J e^^T^ 



log 



_p2 + ]vf2 + 2iE{k)po + 2pk 
-P2 + M2 + 2iE{k)po - 2pk J e^^W ± 1_ 



., (4.11) 



where Im/(po) = [f{po) — /(~Po)]/(2«). We note that the integrals appearing here 
obtain their dominant contributions from the region k ^ T, implying that we may 
replace the logarithms in the unperturbed on shell limit by 

' P"^ + + 2ikpQ + 2pk' 



Im log 



P^ + M^ + 2ikpo - 2pk 



-i log 



Wo+P 
ipo-p 



Im 



E{k) + ipo 

m 



log 



_p2 + ]vf2 + 2iE{k)po + 2pk 
-P2 + M2 + 2iE{k)po - 2pk 



-i log 



ipo+P 
ipo-p, 



(4.12) 
(4.13) 
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as well as set M — > inside the thermal distribution functions. 

Substituting now the above expressions to Eq. ( [4.11|) , we obtain 

VoHP) = ^ log f '-^) r dkk ( + ) (4.14) 

ig'^CpT'^ , fipo+P 
log 



8p \Wo — P , 

from which the value of v'^[P) can be deduced as well. Rotational invariance namely 
implies that 

V (P) = — ^ p = p, (4.15) 

where the ■ P term in the numerator can be read off from the M <^ k limit of 
Eq. (|4lO|) , 

P-yX'^^.P^o = \9'C,T\ (4.16) 
Putting everything together, we obtain to leading order 

V- = ^S^CT' (l - log ( ^) ) p. (4.17) 

Ap^ V 2p \iPo-p// 

It is interesting to note that in the limit under consideration = V~ , implying that 
both bosons and fermions contribute equally to the vector V^{P). 

4.2 Soft masses and momenta 

Let us begin the assembling of the dispersion relation from the case of massless fermions, 
M = 0, and three-momenta of order gT. In this regime, the dispersion relation as solved 
from Eq. ( |4.7| ) clearly becomes (with po i£) 

£-iVo = ±{p+'^\, (4.18) 



P 



from which we further obtain 



m'p +{S±p){m^ log ( I — - ) ± 2pH = 0, (4.19) 



with m? = g'^CpT'^ /2. The two solutions of this equation, corresponding to the + 
and — signs above, are plotted in Fig. ^, where they are compared to the dispersion 
relations of free particles with masses m and 0. 
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Sip) 
m 




^ P 
3.0 



Figure 3: Plots of the dispersion relation of a massless fermion in the SYM theory, <f (p), as a 
function ofp/m with m = ^Cp l2gT . The solid lines correspond to the particle (upper curve) 
and plasmino (lower curve) modes, while the dashed lines show the free massive (upper curve) 
and massless (lower curve) dispersion relations, S-iji) = \/jP^-\-rri? and £{p) = p, respectively. 
In the large p limit, the particle excitation curve approaches the form £{p) = ^/p^ + 2m?. 



The form of the dispersion equation ( [4.19[ ) is equivalent to what has been found for 



quarks in QCD previously [@], except that the effective mass parameter m appearing in 
the result is two times larger in the SYM theory. The two solutions corresponding to 
the minus and plus signs in the above equation are commonly dubbed as the so-called 
particle and plasmino modes, respectively. Of them, only the former one develops into 
a regular particle excitation at larger momenta, while the latter can be identified as a 
collective excitation that is negligible for p ^ gT. This can be verified by evaluating 
the residues of the two poles of the fermion propagator, which can be related to the 
creation probabilities of the respective excitations (see e.g. Ref. [Q). 

The case of fermionic excitations with non-zero, but soft masses M ~ gT can on 
the other hand be studied by numerically solving Eq. ( [4.7|) for various values of the 
mass parameter. The results turn out to be highly similar to the QCD case, which has 
been investigated in Refs. ||, ^ , and we omit the details of the calculations here. The 
main conclusion drawn from them is that with increasing M, the form of the dispersion 
relation of the particle excitation very quickly approaches that of a free particle with a 
shifted mass, 

S = y/p-i + m2, (4.20) 
m = i |m + + 4m2|, (4.21) 
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and already for M = m, this approximation practically coincides with the numerical 
solution of Eq. (|47|) . 



4.3 Hard masses and momenta 

Consider now the case, where either the bare mass M or the external three-momentum 
p of the Af = 2 quark is of order T. In this regime, clearly become of order g"^, so 
the terms in Eq. ( [4.7|) proportional to can be neglected to yield a simpler equation 
for the dispersion relation, 

+ 2P -V + = 0. (4.22) 

Using our earlier result for P ■ from Eq. ( |4.1CI| ) , we see from here that the dispersion 
relation now becomes equal to that of the M = 2 scalars and reads 



S = . /p2 + M2 + 5 Ml (4.23) 



where 



6M^ = 6M. 



2 



1 /""^ „ , 2e^'= 1 k 2e^^^^^ 



7r2 L e^P^ - 1 7r2 



2g'Cp - dkk-^ + - dkk—--:,^^,^ . (4.24) 



E{k) 62/3^^=) - 1 



The behavior of this function has been studied already in Section 



5. Quarks in QCD 

Next, we move on to study fundamental quarks in QCD. To leading order in the 
coupling, the quark self energy can be obtained by evaluating the single graph of 
Fig. 0.C. Continuing with the Feynman gauge, we obtain for this function 

E(P) - 2g C,^^ i2nY i^P + KY + M^R-^ ' ^^'^^ 

which can clearly be written in the form S(P) = [i^ (P) + {if + 2M)a{P)) /2, where 
V~^{P) and a{P) stand for the functions introduced in the previous Section. The 
corresponding dispersion relation is then given by the solution to 

(P(l + a(P)/2) + y+(P)/2)2 + M\l + a{P)f = 0, (5.2) 

which we now proceed to investigate in various limits. The case of soft masses and mo- 
menta has already been worked out in the literature in Refs. ||, H, |^ , but we reproduce 
it below as well for completeness. 
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5.1 Soft masses and momenta 

Let us start from the case of unbroken chiral invariance, M = 0, and three-momenta 
of order gT. In this hmit, the calculation reduces to that presented in the previous 
Section, as we note that we may again altogether neglect the function a{P). The 
dispersion relation is then given by the solutions to the equation 

£ - 2^/2 = ±(p + ^1 , (5.3) 



2p 



from which we again solve 



2m^p+{£±p) jm^ log (f^p^) ±2p^| = 0. (5.4) 

This time the mass parameter fh is, however, smaller by a factor of 2 in comparison 
with the SYM case, 

= g'^CFTyS. (5.5) 

The interpretation of the two modes corresponding to the + and — signs above is 
equivalent to the SYM case, and the behavior of S{p) is trivial to infer from Fig. |^. 
The case of non-zero but soft masses has been studied in Ref. 0, with the result that 
the collective plasmino mode quickly becomes negligible when the rest mass parameter 
exceeds the value in. 

5.2 Hard masses and momenta 

Next, let again either the bare mass or the external three- momentum be of order T, 
which to the best of our knowledge has not been studied in the literature before. In this 
limit, the functions a{P) and V~^{P) are both regular and of order g'^ in the unperturbed 
on shell limit, so we have from Eq. (|5.2|) the dispersion equation 



= E{pY+ P-V+{P) + M^a{P) . (5.6) 

L J p2=-M^ 

Using the results presented in Section ^, we immediately obtain for the first integral 



a result obviously independent of the three-momentum. Similarly, we get for the func- 
tion a{P) 

g'Cpj 1 f P' + + 2zkpo + 2pk \ 1 

"^^) = ^r^yo '^''p''^[ p^ + M^ + 2^kpo-2pk )^^ ^'-^^ 

^ r ,1 ^ , f-P^ + M^ + 2iE{k)po + 2pk 
+ Re dk log ' 







pE{k) V + M2 + 2iE{k)po - 2pk J e^^^C^) 
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e{pf-E{pf 



0.24 
0.22 





Figure 4: The thermal correction to the non-interacting dispersion relation of a quark in 
QCD, S{p)'^ — E{p)'^, as a function of p/T for various values of M/T. Counting from top 
down on the right side of the plot, the curves correspond to M/T = 0, 1/2, 1, 2, oo. 



in which the first logarithm vanishes upon taking its real part and then going to the 
on shell limit. From the second logarithm, we on the other hand get 



where we have used the result of Eq. (|3.4| ) . 

First consider the case of M = but p ~ T ^ gT. In this limit, M'^a{P) naturally 
vanishes, while P ■ V^{P) = ^g'^CpT'^, so the dispersion relation reads 



S{pY=p^ + m^, (5.10) 

with = jg'^CpT'^- Similarly, in the case with M ^ cl{.P) as well as the second 
term in Eq. ( |5.7| ) become exponentially small, so the dispersion relation becomes 

^(p)2 = p2 ^ + (5.11) 



with 

= -g^CpT^ (5.12) 
6 

Between these two limits, we must resort to a numerical evaluation of the functions 
P ■ V^{P) and a(P), and obtain the dispersion relation £{j)) displayed in Fig. |[ 
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Finally, we note for completeness that if M is of order T, but p ~ gT, we may 
obviously expand Eq. ( p.9|) in the limit of small p/k. Using Eq. (|5.6| ), this leads to the 
dispersion relation 

S{p)^ = p2 ^ + 6M^, (5.13) 
SM^ = ^^—^ / dk < — + — } , (5.14) 

where SM"^ corresponds to the p ^ limit of the functions displayed in Fig. ^ 



6. The non-relativistic heavy particle hmit 

As we have seen above, the leading 0{g'^) thermal corrections to the dispersion relations 
of heavy fundamental quarks and scalars are suppressed by an inverse power of the bare 
mass and consequently vanish as M ^ oo. At the same time, on physical grounds one 
would expect the thermal corrections coming from soft loop momenta of order gT 
to survive. To see this, consider the interaction of a heavy charged particle with an 
electromagnetic field, and in particular, how thermal effects modify the long distance 
behavior of the fields sourced by the particle. Debye screening cuts the longitudinal 
electric fields off at a distance of order the inverse Debye mass rriu, which can easily be 
seen to give the energy stored in the electric fields a temperature dependent contribution 
of order g'^ruD- The factor g'^ here comes from the fact that the electric fields sourced 
by the quark are 0{g), so the electromagnetic energy density (which depends on the 
square of the electric fields) is O^g"^), while the factor m^) is a direct consequence of 
the finite range l/mo of the static fields.^ We therefore see that at cubic order in the 
coupling, one should find thermal corrections to the heavy particles' energy (and thus 
their dispersion relation), which are independent of the heavy mass and in fact become 
the dominant contributions for bare masses of order T/g and larger. 

To put the above statements on a more quantitative footing, consider now the 
QCD type diagram of Fig. Q.c, which stands for a generic example of the one loop 
graphs contributing to the heavy fermion/scalar self energies in both QCD and SYM. 
The IR sensitive region of momentum space that might upon a suitable resummation 
lead to the suspected 0{g^) behavior of the quantity is one where the gluon line is 
soft, i.e. qo ^ q ^ gT. To inspect which types of loop corrections the graph may 
be sensitive to, we add an extra gluon line to dress either the light field propagator, 
the heavy field propagator, or the vertex, as displayed in Fig. ^ and investigate, what 

^Due to the Debye screening, the energy stored in the fields becomes proportional to an integral of 
the form f d^p ^} . 
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b) ^ c) d) ^ 

Figure 5: Four different Feynman diagrams depicting possible two-loop corrections to the 
heavy quark self energy in QCD. As is explained in the text, only the graph a is observed to 
lead to an 0{g'^) contribution to the quantity. 

power of g these corrections are suppressed by. It is easy to convince oneself that the 
qualitative results obtained here will also apply to the other diagrams of Fig. |l].a-b 
with the same topology, where either a gluon or a light scalar line is dressing the heavy 
quark or scalar. A brief inspection on the other hand shows that the graphs involving 
either a light fermion propagator or a four boson coupling cannot contribute to the self 
energy at the order we are interested in. 

For the graph of Fig. |^.a, which is a generic example of diagrams containing loop 
corrections to the gluon propagator in Fig. |l|.c, we observe that the maximal contri- 
butions originate, when the four-momentum of the added propagator is hard, K T. 
Comparing this to the original one-loop graph, we observe that it obtains an extra 
factor of from the two new vertices, as well as a factor 1/(7^ from the induced soft 
gluon propagator with momentum ~ gT (recall that we assumed the momentum flow- 
ing along the gluon line of the original graph to be soft). We therefore conclude that 
the addition of the extra gluon line contributes at relative order g^ and that loop cor- 
rections to the soft gluon line of the graph of Fig. |I|.c must thus be resummed. This 
resummation may be carried out using the hard thermal loop (HTL) approximation 
for the gluon self energy, as the momentum flowing along the gauge boson line of the 
original one- loop diagram is soft, while the momentum flowing in the loops of the self 
energy is hard. 

For the graphs b-d, it is straightforward to see that similar order-of-magnitude 
estimates, both for hard and soft loop momenta, always lead to the suppression of the 
result by at least one relative power of g. This behavior can be seen to carry on to 
further perturbative orders as well, thus confirming that it will not be necessary to 
carry out resummations amounting to dressing the heavy field propagator or the vertex 
function in the original diagram. We thereby see that in order to obtain the 0{g^) 
contributions to the self energy of a heavy quark in QCD, we only need to consider 
the graph of Fig. ^.c, where the gluon line has been dressed with the HTL self energy. 
This diagram, as well as the corresponding ones in Figs. p.a and b, needed in the SYM 



computation, will be evaluated below in Section |6.2| . In the graphs containing a light 
scalar exchange, the HTL self energy reduces to a scalar thermal mass squared, m^. 



-17- 



a) 




,-0-, b) 




Figure 6: The diagrams contributing to the self energy of an infinitely massive fundamental 
a) scalar and b) fermion in the = 4 SYM theory, and c) a quark in QCD. The intermediate 
light propagators are dressed with the HTL self energy function, which for the light scalars 
implies the use of the scalar thermal mass m^. In the limit of a large quark mass AI > T/g, 
the graphs are dominated by soft light quark momentum go ~ '7 qT- 

Before proceeding to the resummation calculation, we note that the above reason- 
ing also implies that it is possible to determine the thermal corrections to the dispersion 
relation of a heavy particle using classical field theory. As the momentum flowing along 
the dressed boson lines in Fig. ^ is parametrically smaller than T, thermal fluctuations 
of the light fields sourced by the heavy particle should be negligible and the dynam- 
ics of these fields governed by classical field theory. The necessity to dress the light 
propagators with the HTL self energy, however, implies that the fields will experience 
thermal effects — the classical equations of motion and the classical Hamiltonian will 
depend on thermal susceptibilities corresponding to the HTL self energy function. For 
a gauge field, the classical equations of motion will be those of macroscopic electrody- 
namics with temperature dependent permittivity and permeability tensors, which can 
be straightforwardly obtained from the HTL self energy function. For a scalar field, the 
equation of motion will on the other hand simply be the Klein-Gordon equation with a 
temperature dependent screening mass, which is nothing but the zero momentum limit 
of the light scalar field self energy function. 

We will now proceed to consider the classical calculation outlined above, with the 
objective being the derivation of the correct large-M limit of the dispersion relation of 
the fundamental particle. After this, we will in Section verify the result by perform- 
ing the corresponding diagrammatic computation, explicitly evaluating the graphs of 
Fig. ^. As we are here considering heavy particles with typical momenta p ~ T <^ M, 
we will in both cases work in the non-relativistic limit and thus restrict ourselves to 
determining the thermal shifts to the rest and kinetic masses of the excitation. 

6.1 Classical field theory 

With the above arguments in mind, consider now a large volume V, which encloses a 
heavy charged particle moving along some trajectory. The energy contained in V will 
depend on how fast the particle is moving. As we are taking M — oo, it is convenient 
to consider the dispersion relation as a function of the particle's velocity v rather than 
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its momentum. Classically, the energy contained in V is 



M 

£{v) = _ + A£:fie,d(^;), (6.1) 



V 



where A£^ficid is the shift in the energy stored in the light scalar and/or electromagnetic 
fields sourced by the charge due to screening efi^ects.^ By rotational invariance, the 
field energy must have the small v expansion 

A£:fie,d(^;) = <^M,est + ^SM^^y + 0{v^). (6.2) 

This equation, which is manifestly finite in the large M limit, defines the shift in the 
rest and kinetic mass of the particle. 

In what follows, we will compute the contributions to the rest and kinetic masses 
of the heavy fundamental particle from the light scalar and electromagnetic sectors 
separately, and then assemble the results for quarks in QCD as well as quarks and 
scalars in SYM. We will henceforth refer to the heavy particle as a quark, but this is no 
restriction, as the classical calculation is not sensitive to the spin of the heavy particle. 

6.1.1 Scalar field contributions 

Let a; quark (^) now denote the trajectory of a heavy quark with mass M. In theories 
such as A/" = 4 SYM, the quark's coupling to the light scalar fields is dictated by an 
effective Lagrangian of the type 



^classical 



= -M^l - x^„3.k + / d'x \l (0^ - (V<^)^ - my) 



P<t> 



(6.3) 



where 



P = eeff-^l-i^uark ^^{^ - iKquark(i)) (6.4) 



is the charge density of the quark, which must be a Lorentz scalar. The constant egg 
denotes the effective coupling of the fundamental particles to the light scalar fields, and 
through matching to the full theory^ is easily seen to read el^ — CFg^- The parameter 
rUs on the other hand corresponds to the thermal screening mass of the light scalars, 
and is thus a function of the temperature. Its value is given by 

m^ = hmn,(Q), (6.5) 



^The renormalized mass M already contains the energy of the fields in the vacuum. 

^An easy way to see that this is the correct value of eefi is to compute the energy density stored 
in scalar fields at distances d <C 1/T from the quark and match the resulting expression to the full 
theory result. 
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where Ils{Q) is the hght scalar self energy function and the limit of vanishing four- 
momentum is independent of the order of taking pq and p to zero. 

From the Lagrangian of Eq. (|6.3|) , we easily obtain the classical Hamiltonian 



H, 



M 



classical 



scalar ; 



(6.6) 



X 



quark 



where 



scalar 



+ {Vcpf + m 



X 



quark 



(6.7) 



is the energy stored in the scalar field. This leads to a sourced Klein-Gordon equation 
of motion for the field, 



(-92 + m^) = p, 



(6.8) 



which, taking the quark's velocity v = iquark to be constant and introducing a spatial 
Fourier transform, has the solution 



0(t, k) = eeffVl 



-iu!o{k)t 



(6.9) 



with Uo{k) = V ■ k. Substituting now this into Eq. (|6.7|) , we see that we can finally 
write the classical field energy in the simple form 



^scalar „ 



d^k 2ujI - v'^ {ujI + k'^ + m]) 



i2n) 



-UJn 



(6.10) 



which we however immediately observe to be UV divergent. 

To cure the divergence, we note that it is independent of the screening mass 
and hence the temperature T, and corresponds to the singular ~ l/x^ energy density 
of a localized point particle. We may thus simply subtract from the field energy the 
corresponding expression with nis = 0, 



AS, 



scalar 



scalar 



scalar 



ms=0 ' 



(6.11) 



which is allowed as we are only interested in the shift in the field energy due to the 
presence of screening. The resulting three-dimensional momentum integral is both 
finite and straightforward to evaluate, and leads to the simple result 



AS, 



scalar 



-eff 



1 



(6.12) 
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Taking the small v limit here and comparing to Eq. (|6.2|) , we see that the shift in both 
the rest and kinetic masses of the heavy quark due to its interaction with the light 
scalar fields is equal to 

8M = ^S^. (6.13) 
Stt 

This equality is a consequence of the Lorentz invariant nature of the influence of Iliight 
on the classical equations of motion, as in the long wavelength limit it reduces to the 
Lorentz preserving screening mass mg- As we will discuss below, this is however not 
the case for electromagnetic fields. 

6.1.2 Electromagnetic field contributions 

The derivation of the energy contained in electromagnetic fields is slightly more in- 
volved than it was for the case of scalar fields. The reason is that even in the limit 
of a soft external momentum, the gluon polarization tensor has nontrivial momentum 
dependence, which results in complex permittivity and permeability tensors, and cor- 
respondingly non-trivial dispersion and damping of the electromagnetic fields. The 
presence of dispersion and damping makes the construction of a Hamiltonian from 
the Lagrangian formulation of the electromagnetic problem impossible. However, as 
we elaborate on below, using the classical equations of motion and the expected form 
of the energy conservation equation in the presence of dissipation, one can extract a 
meaningful definition of the electromagnetic energy of a moving quark. ^ 

We begin the calculation by noting that in the limit of weak fields, the non-linear 
non-Abelian interaction terms appearing in the classical Yang-Mills equations of motion 
may be neglected. This leaves as the linearized equations of motion of the electric and 
magnetic fields the macroscopic Maxwell's equations in a medium, 

V-D = p, V-B = 0, (6.14a) 

dB dD 
\/xE = - — , VxJ/ = J + — , (6.14b) 

where the macroscopic fields D and H are defined in Fourier space [ko = u here) via 

D{iu, k) = e{uj, k)E{u, k), (6.15) 
B{iu, k) = fi{uj, k)H{uj, k), (6.16) 

and e and /x are the permittivity and permeability tensors describing the medium. We 
note that even in the limit of small momentum, e{uj,k) and fj,{uj,k) still have nontrivial 



^The argument we use below is a straightforward generalization of that found in Refs. |lT| . 
Similar issues have also been discussed in Ref. lOl. 
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momentum dependence, as they become functions oiuj/k. The source = {p,J) is 
on the other hand simply the current of a massive localized quark, 



p = CcsS {X - ajquarkW), (6.17) 

~dt 



J = e,s^^^6-\x-x^,^,^{t)). (6.18) 



Consider now a large volume V, which encloses the source and whose size is 
larger than any characteristic length scale associated with the medium. In this volume, 
the energy conservation equation reads 

^^rf^XMEM = - ^rf'x{V-5+ J-Sj-Qheat, (6.19) 

where uem is the electromagnetic energy density, S = E x H the Poynting vector, 
and Qheat the energy lost to heat per unit time. This equation simply states that the 
rate at which the electromagnetic energy contained in V decreases is given by the rate 
at which electromagnetic energy flows out of V plus the work done by the source and 
the energy lost to heat. Using Maxwell's equations, we can furthermore express the 
divergence of the Poynting vector in the form 



V.S=-£;.f -H.f -J.£. (6.20) 



flnally giving us 



In a medium where there is no dispersion or dissipation [i.e. where e and ^t are 
constant), the time derivatives appearing on the right hand side of Eq. ( |6.21| ) can be 



moved out of the integral. In this case, we furthermore have Qheat = 0, and thus 
uem = \ {E ■ D + B ■ H) . If e and ^ are not constant, as is the case in the theories 
we are considering, the time derivative, however, cannot be factored out of the inte- 
grand, and more care must be taken in order to obtain the correct expression for the 
electromagnetic energy. 

To this end, let us now introduce a spacetime Fourier transform and assume that 
the electromagnetic flelds may be written in the form 

E{uj,k) = E{uj,k)f{uj -ujQ{k)), (6.22) 
H{uj,k) = H{u,k)f{u-ua{k)), (6.23) 
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where the real valued function /(x), normalized such that dx f(x) = 1, has support 
only on a narrow interval around x = 0, and E{uj, k) and H{uj, k) are slowly varying 
functions around u = ujQ^k).^ Under this assumption, we obtain 



duj duj' d^k 



(6.24) 



+ 



27r 27r (27r)3 
dui dui' d^k 
2^ 27r (27r)3 



-2cj£;t(c^', fc)e(cj, fc)^;(cu, fc)} /(cj' - oJo{k))f{uj - cJo(fc))e^^^'"^^* 
-iujW{uj', k)^{uj, k)H{uj, k)} f{uj' - LOo{k))f{uj - cJo(fc))e^("'-")*, 



where the frequency integrations will be dominated by uj 



~ /.i' ~ 



UJ 



ujQ{k). Expanding 



now the slowly varying part of the integrand in powers of {u — uj') — which can be 
represented as a time derivative — and integrating over the frequency variables, we 
find that at leading order the result reads 



d^x I E 



V 



where we have defined 



dD 



+ H 



dB 



d& 



EM 



dt 



+ Q 



heat ) 



(6.25) 



^EM = 



d^k 



E\uj,k)d^[ k)]Eiuj,k) + H\uj,k)d^[ ujfi{uj, k)] H{uj, k) 

(6.26) 



87r2 J (27r)3 

+ ujE^{lj, k)e{uj, k) d^E{uj, k) - uod^E^uo, k)e{uj, k)E{uj, k) 



+ ujH\uj, k)fi{uj, k) d^H{uj, k) — ujd^H\uj, k)fi{io, k)H{uj, k) 



uj=uJo{k) 



and 



Qheat — 



An' 



d^k 
(2^ 



iuoE^uo, k) e(a;, k)E{LJ, k) 
i(jjH\uj, k)fi{uj, k)H{uj, k) 



(6.27) 



uj=uJo{k) 



Comparing the above expressions to the anticipated form of the energy conser- 
vation equation ( p.21| ), we identify the quantity Sem with the electromagnetic field 
energy contained in V. In the limit, where e and fi are constant, only the first 
two terms in Eq. (|6.26|) contribute to the integration, and as expected, we obtain 



^Such an ansatz is certainly valid for fields sourced by a quark moving at a constant velocity in 
which case they become proportional to S{uj — v ■ k). 
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^EM = I / {E ■ D + B ■ H). The quantity Qheat, which only depends on the imagi- 
nary parts of e and /x, is on the other hand clearly the energy lost to heat in the volume 
per unit time. 

To obtain the shifts in the rest and kinetic masses of the quark, we now again 
specialize to the case, where its motion non-relativistic. In this limit, the magnetic 
field and the transverse component of the electric field created by the quark are 0{v) 
already at T = 0, while thermal contributions to /x and the transverse component of 
e are 0{uP) in the small momentum limit. We therefore conclude that any thermal 
effects in the electromagnetic field energy coming from magnetic and transverse electric 
fields are 0{v^), and we may concentrate on the contribution to the field energy coming 
from longitudinal electric fields. 

Next, we solve for the longitudinal component of the electric field from Maxwell's 
equations by first introducing a spacetime Fourier transform. This leads us to 

El{u, k) = 27reeff f ,M uj - uj^{k)), (6.28) 

where again Uoik) = v ■ k, and is the longitudinal part of the permittivity tensor, 
which can be given in terms of the gluon polarization tensor as 



eAuJ^k) = 1 + 



fc2 



In the HTL limit, we on the other hand obtain for IIoo 

' k + oj 



H 13 



m 



D 



1-^ 
2k 



log 



k — uj 



ITX 



(6.29) 



(6.30) 



where m^) is the Debye mass of the theory. 

Putting finally Eq. ( |6.26| ) and Eqs. ( |6.28| )-( |5.30| ) together and expanding the re- 
sulting expression in powers of f , we find that to second order in the velocity, the 
electromagnetic field energy is given by 



_fk_ 

(2^ 



{Ak^ml - (vr^ - 4) m%) 



2 ^3 
d) 



V 



kf 



(6.31) 



As in the case of scalar fields, the above integral contains a UV divergence independent 
of the Debye mass, which we remove by considering the difference 



-EM 



ni£)=0' 



For this quantity, we straightforwardly find 



(6.32) 



(6.33) 
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where 



5M,est = (6.34a) 

2 ™ / ^2 



As previously advertised, the contributions to the rest and kinetic masses of the 
heavy quark coming from the electromagnetic field are not equal. From the above 
results, we see that the shift in the kinetic mass is smaller than the shift in the rest 
mass, and note that the discrepancy can be traced back to the imaginary part of 
the permittivity tensor, which is itself related to the presence of dissipation in time 
dependent electromagnetic fields. Evidently, the presence of dissipation decreases the 
energy stored in time-dependent electromagnetic fields. 

6.1.3 Assembling the results 

With Eqs. ( |6.13D and ( |6.34| ) at hand, we are now finally ready to assemble our results 



TT 



for the thermal corrections to the dispersion relations of heavy quarks and scalars in 
SYM and heavy quarks in QCD. In the SYM theory, there are six light scalar fields for 
each gauge boson. From the Lagrangian of Eq. (|2.4| ), we however see that the M = 2 
fields couple to only one of these, bringing the shifts in the rest and kinetic masses to 

dMrest = 3 \ z ' (6.35a) 

OTT OTT 



Here, we have expressed Ceg in terms of the group theory factor Cp and the full theory 
coupling g, while the thermal masses appearing in the results are given by 

ml = 2N^g^T^ = 2ml (6-36) 

For QCD, where the gauge fields are the only light bosonic fields, the shifts in the 
rest and kinetic masses on the other hand read 

5M.3t = (6.37a) 



where the Debye mass has the value 

i2t^2 



ml = ^-^ (n^ + ^ ]. (6.38) 
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As in the SYM case, we note that this result represents the leading thermal correction 
to the dispersion relation of a heavy particle for bare masses M > T/g. For smaller 
mass scales, the leading order result can be read off from Eqs. ([4.24|) and ( p.llj ). 

6.2 Quantum field theory 

To complete the discussion from the beginning of this Section and to verify the result 
of the classical calculation, we will next compute the leading order thermal correction 
to the dispersion relation of infinitely heavy fundamental particles using the techniques 
of perturbative quantum field theory. We will do this only for an A/" = 2 test quark 
or scalar, for which we are lead to consider the graphs of Fig. |l|.a-b with one gluon 
or light scalar line dressing the massive quark or scalar. The QCD case can later be 
trivially obtained from the M = 2 quark result by simply setting the scalar thermal 
mass to zero. 

If we take the bare mass of the heavy quark to be infinitely large, we obviously 
must first take this limit and only then perform an expansion in powers of the coupling 
constant. A simple power counting exercise shows that for quarks and scalars with 
masses M > T/(yf, the dominant contribution to the graphs originates from the region, 
where the momentum flowing along the light boson line is soft, qo ^ q ^ gT. As 
explained above, this implies that we must use a gluon propagator dressed with the 
hard thermal loop (HTL) self energy and a light scalar propagator containing the scalar 
thermal mass m^, leading us to the graphs of Fig. ^a-b. 

Heavy fermion dispersion relation 

Let us first consider the heavy M = 2 fermion dispersion relation in the SYM theory, 
and thus the two graphs of Fig. ^.b. Denoting the value of the first diagram by Sfps(P) 
and the second by S^gg(P), we obtain for the resummed heavy quark self energy 

Sre.(P) =SL(P) + S^es(^), (6.39) 

where 

T.UP) ^ g^C.i'^j^^^^^'i^D-'iQ), (6.40) 

and D'^" denotes the resummed gluon propagator. 

Performing now some straightforward algebra and using the fact that to leading 
order, the quark satisfies the Dirac equation, we can identify the leading large-M term 
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of the gluonic part of the self energy as 

In Appendix A, wc show that to leading order in a large-M expansion, this function 
can be evaluated to yield 

E,es(P) = 2g^CFE{p)hjo, (6.43) 
where the real part of the integral is given by 

For the scalar contribution to the self energy, we on the other hand obtain 

To 



ml 



/$((P-Q)2 + M2)g2(Q2 + ^2) 



8. M^ + ^ilFj- 

Here, we have first regulated the UV divergence of the integral by subtracting from 
it the corresponding unresummed expression, and then used the fact that its leading 
large-M piece originates from its simple T — hmit. 

Using the above results, we may now solve for the dispersion relation of the heavy 
quark from the poles of the dressed fermion propagator + M + Eres(-P)), which 

leads us to the equation 



Att Att M 



1 



3M2 V 64 



■. (6.46) 



Solving S{p) from here and expanding the result in the limit of small velocity, we finally 
obtain as the leading order rest and kinetic masses in the large M limit'' 

M.e. = M{l-/C.!^^}, (6.47) 

Mkin = m|i - q^Cf^^ g^Cp-^ \. (6.48) 

\ ^ 247rM V 16 y ^ SttM j ^ ' 

These results arc seen to be fully consistent with those of the classical computation. It 
is also worth noting that the imaginary part of II does not contribute to the result at 
this order yet. 

^Note that here the energy is expressed in terms of the canonical momentum of the quark, and thus 
the rest and kinetic masses are defined through the effective Hamiltonian, H = Mrest +P^/(2Mkin). 
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Heavy scalar dispersion relation 

For the case of a heavy M = 2 scalar, the diagrams we need to consider are depicted 
in Fig. p.a. In the large mass limit, these graphs lead to the self energy reading 



and ultimately to a dispersion equation of the form 



Att M 



1 - 



2v' 



1 - 



TT 

64 



(6.49) 



(6.50) 



Comparing this to Eq. ( |6.46|) , we observe perfect agreement with the case of a heavy 
M = 2 fermion. This was of course to be expected due to the unbroken supersymmetry 
of the theory in the limit under consideration. 



7. Discussion and conclusions 

In this paper, we have considered thermal corrections to the dispersion relations of 
fundamental particles both in an A/" = 4 Super Yang- Mills plasma coupled to an A/" = 2 
hypermultiplet and in the quark gluon plasma (QGP). We have performed the calcu- 
lations to the leading order in a weak coupling expansion, our main motivation having 
been the study of the nature of the excitations at various values of the bare (zero 
temperature) masses and three-momenta and a comparison of the results between the 
two theories. The fact that we have left out higher order corrections from our studies 
implies that our results are applicable only in the limit of very weak couplings (high 
temperatures in the case of QCD), and that we are not able to study the decay widths 
of the excitations. 

Our results for the dispersion relation in different kinematical regions of the two 
theories can be read off from Eq. ( p.5| ) for the Af = 2 scalars, Eqs. ( [4.19| ) and ( |4.24| ) for 



the M = 2 quarks, and Eqs. ( p. 4] ) and ( p. 61 ) for quarks in QCD. They are summarized 
in Table 1 displayed at the end of the paper. Inspecting the results and comparing 
them to each other, some interesting patterns can now be observed. 

If the bare masses of the particles are small, the dispersion relation of the M = 
2 quarks behaves in a way closely analogous to the dispersion relation of quarks in 
the QGP. For vanishing bare masses and soft three-momenta, the quantity is in both 
theories obtained from an equation of the form 

2m^p + (^ ± p) jm^ log (j^^ ± V| = 0, (7.1) 
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where only the expression for the soft mass parameter m ~ gT differs between the 
theories. This equation exhibits two solutions, a particle and a collective plasmino 
excitation, which are analyzed in some detail in Section 4. In the limit of larger three- 
momenta, p ~ T, the particle solution reduces to a simple form 

S^=p' + ml, (7.2) 

where the asymptotic mass rrioo is related to m above such that = 2m^. 

The above results for fermions are in stark contrast with the dispersion relation of 
the M = 2 scalars. For the latter, the dispersion relation is for all values of the bare 
mass M and three-momentum p given by 

£ = ^v" + M2 + 5M|, (7.3) 

where (5M| is independent of the momentum p but depends on the bare mass M as 
shown in Fig. 2. We thus see that for soft bare masses and three-momenta, the behavior 
of the dispersion relations in the SYM theory does not exhibit supersymmetry: The 
fermions and scalars within the same SUSY multiplet behave differently. This can be 
attributed to the presence of a non-zero temperature in the system, which is known to 
explicitly violate SUSY. 

Approaching bare masses of order T, the dispersion relation of the quarks can in 
both theories be written in the form 



£= ^p^ + M^ + 5M^. (7.4) 

Here, the mass correction is equal to the scalar mass correction SM^ for the M = 2 
quarks, while for quarks in QCD, it becomes a slowly varying function of p that we 
displayed in Fig. 4. We thus observe that the dispersion relations of the M = 2 quarks 
and scalars become equal in this limit, and that supersymmetry is restored at least in 
this sense. The same is true, if the bare mass is small but the three-momentum hard, 
since SMg approaches moo in this limit. 



The above observations are in agreement with the conclusions drawn in Ref. [15|, 
and can be easily understood in the limit, where either M or p becomes the dominant 
scale in the problem. Although the heat bath explicitly breaks SUSY, the breaking 
of the symmetry is not manifest in the propagation of the particles unless they are 
sufficiently light and soft, as heavy Af = 2 propagators act in the diagrams essentially 
as zero-temperature ones. For fundamental particle masses M >T, the leading thermal 
contributions to the self energies in fact come from the screening of the Af = 4 fields, 
to which the fields in the M = 2 hypermultiplets couple to identically, with couplings 
independent of T. What we find to be somewhat surprising is that in our calculations 
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this effect is seen to take place as soon as M or p is parametrically larger than gT, 
and it is only strictly in the limit of soft momenta and bare masses that we see any 
difference between the M = 2 fermions and scalars. 

Moving on to yet larger masses, a particularly interesting special case of our results 
is encountered, when the bare mass of the fundamental fields becomes parametrically 
larger than the temperature. In this limit, we may compare our results to the findings 
of Ref . [Q , where the authors considered the non-relativistic dispersion relation of heavy 
(M ^ \/~XT = y/T^gT) M = 2 quarks and scalars in the large A*";, limit of a strongly 
coupled SYM plasma. According to their results, the excitations follow a dispersion 
relation of the form 



p 



(7.5) 



2Mkin(T) 

where both Mrest and Mkm can in the limit of a large rest mass be given in the form 



m ( 



m 



(7.6) 



(7.7) 



and the difference between Mrest and Mkin is suppressed by at least one inverse power 
of M and is displayed numerically in Fig. 1 of Ref. [Q. 

This behavior is to be contrasted with our results from Section |^, which state that 
for heavy non-relativistic M = 2 quarks and scalars immersed in a weakly coupled SYM 
plasma, the shifts in the rest and kinetic masses are given by Eq. ( 6.35| ). Going to the 
large limit and defining 



m = 



Sir 



we obtain from Eq. ( |6.35 ) 
Mrest = m|i- 



M 



kin 



M< 1 




(7.8) 

(7.9) 
(7.10) 



We observe that the sign of the shift in the rest mass is the same both in the weakly 
and strongly coupled limits: In both cases thermal effects lower the value of the mass. 
However, the signs of the shifts in the kinetic mass do not agree. In the weakly coupled 
limit, the contribution to the quantity from gauge fields is negative, but smaller in 
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magnitude than the effect of the scalar fields, which gives the shift a positive overall 
sign. Moreover, in the weakly coupled limit the difference Mj-est — M]^i^ is constant and 
non-zero as M ^ cxd, in contrast to the result in the strongly coupled limit. 

Some interesting observations can now be made. Although the shifts in the ef- 
fective rest and kinetic masses are not equal in the weakly coupled limit, the scalar 
contributions to them are, implying that the difference in the masses arises from the 
dynamics of the gauge fields alone. At leading order in a weak coupling expansion, a 
physical difference between the scalar and gauge field dynamics is that the former is 
not dissipative. Indeed, it is easily seen that a non-zero difference between the kinetic 
and rest masses requires a momentum dependent permittivity tensor. Moreover, by the 
Kramers-Kronig relations, any momentum dependence in the real part of the permittiv- 
ity tensor implies that its imaginary part must be non-zero. At least at weak coupling, 
we therefore see that there is a connection between having a non-vanishing difference 
between the rest and kinetic masses and the dissipation of gauge fields. While this 
argument need not hold at larger coupling, it is interesting to contrast our observation 
with the fact that the mass difference is suppressed by 1/M in strongly coupled SYM 
theory. We leave the exploration of this subject for future work. 
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A. Evaluation of Sf^^ 

In this Appendix, we evaluate the large-M limit of the sum- integral defined by Eq. ( |6.42| ). 
It is needed in determining the leading order thermal correction to the heavy quark 
dispersion relation due its interaction with a dressed gluon, and is relevant in both SYM 
and QCD. We will perform the calculation in an expansion in powers of the velocity v 
up to and including order f^, as in Section 6 we are ultimately only interested in the 
thermal shifts to the rest and kinetic masses of the quark. 

Taking advantage of the gauge invariance of the quantity under consideration, we 
choose to work in the Coulomb gauge, in which the dressed Euclidean space gluon 
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propagator reads 



+ 



g2 + noo ' q' + ql + Ii^' 
with the components of the HTL polarization tensor having the forms 



(A.l) 



noo(Q) 
nT(Q) 



D 



-m 



1 + ^ln^^ 

2q iqo + q 



D 



(A 



^ mjqp + q^) iqo - q 



4g3 



^qo + q 



(A.2) 
(A.3) 



and rriD denoting the Debye mass of the theory. We also choose the branch cut of the 
logarithm to be slightly tilted off from the negative real axis, so that these expressions 
have a branch cut going from go = —"iq to go = "iq via the left side of the complex plane 
(see Fig. |^). This is required in order to keep the polarization tensor an even function 
of go when evaluated at real Matsubara frequencies, while maintaining the contribution 
of the small- 1 go I limit well-defined. 

Having the vector p point in the z direction of our coordinate system, it is easy to 
see that the above produces 



■ Po7o + 



J-oo 



g2 + g2 + Ht 



P^li . (A.4) 



We regulate the UV divergences appearing here by subtracting from the integrand the 
corresponding unresummed expression, which in addition helps us to isolate the leading 
large-M limit of the graph we are interested in.^ This gives us 



-2ig'^CF 



1 



q" + Hoo 



Polo 



+ 



ql + q^ + Ht ql + q^ , 
= -2z5(^Cf< hPolo + hPilt 



Pili 



(A.5) 



which defines the longitudinal and transverse integrals, II and It, that we now set out 
to evaluate. As the imaginary parts of the integrals only affect the dispersion relation at 
higher orders in g than what we are pursuing in this work, below we will only consider 
their real parts. 
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Figure 7: Illustration of the contour C on the complex qq plane. The red crosses denote the 
locations of the bosonic Matsubara modes at go = 2n7rr, the blue wiggly line the branch cut 
of the integrand going from qq = —iq to qq = iq, and the green curve our integration contour. 
The contour C encircles all the poles of the cotangent function, yet avoiding crossing the 
branch cut. 

A.l The longitudinal integral 

Let us first look at the longitudinal case. The Matsubara sum appearing in the defini- 
tion of the integral can be transformed into a contour integral using standard tricks, but 
recalling that there is now a branch cut in the definition of the polarization tensor that 
the integration path should avoid. Choosing thus a somewhat non-standard contour C 



Here, we have first neglected the part of the cotangent function that vanishes in the 

^It is easy to confirm that for the unresummed integral, the finite temperature contribution is sup- 
pressed by two inverse powers of the heavy mass, while the T = piece is part of the zero temperature 
quantum corrections we are not interested in. We are thus justified to make the subtraction. 



according to Fig. ^ and defining g{qo) 




^ ), we obtain 




(A.6) 
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T ^ limit^, then denoted the point at which the branch cut intersects the real axis 
by — c, and finally used the fact that the function g is analytic on the Euclidean real 
axis to proceed to the e ^ limit. In the last form, we are integrating the function 
g through the branch cut, which is obviously the correct prescription in this case. It 
will be seen later that the location of the branch cut — as long as it doesn't lie on the 
imaginary axis — does not affect the final result of the integral. 
Writing now the integral in its full form, we have 



dS 1 r 1 1 

(Po - qo? + {p- qf + M^W + noo(go, q) ~7 
d^q 1 noo(go,g) 



(27r)4 g2 + ^2 _ 2iqo^/WTf - 2q ■ p q\q^ + noo(go, q)) 



(A.7) 



where in the latter form we have already taken the on shell limit of the heavy particle, 
Pq = iE{p). This implies that in order to correctly perform the analytic continuation in 
the integral, we must slightly deform the go integration contour, so that it goes between 
the two poles of the fermionic propagator, as it did for real pq. 

To further analyze the integral of Eq. ([A.7|) , we note that if either qo or q was 
parametrically larger than rriD, {e.g. of order M), then it would be easy to see that 
the integral would scale as (m£)/M)^+", with a > 0, i.e. be subleading with respect 
to the contributions we expect to find. We will thus rescale the integration variables 
according to 

(A.8) 

q = m£)fc, (A. 9) 

go = rriDkz, (A. 10) 

noo(go/g) = mlf{z) (A.ll) 
p = Mv, (A. 12) 

in terms of which we have 

1,=-^ f ^ /^l I fA±_ (MS) 

Jc'2n J {271^ k x{l + z^)k^-2izkVTT^-2v-k + f{z)' ^ ' ' 

with C reminding us of the deformation of the go integration path. Note that in our 
notation throughout this Appendix, v just a shorthand for p/M, and not the correct 
relativistic velocity of the heavy particle. This is not a problem, as we are merely 
interested in the two first terms of the small-w expansion of the dispersion relation. 



'On dimensional grounds, this is again suppressed by at least two powers of M . 
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Figure 8: Illustration of the contour C. The red crosses denote here the two imaginary 
poles of the fermion propagator, but otherwise the notation is as in Fig. |^. 

Next, we point out that out of the two poles of the fermion propagator 

z = z± = ^ (yTTv^ ± v^l + (v-xk)^^ , (A.14) 

both are imaginary, one being large i) and always on the upper half plane, while the 
other can reside on either side of the real axis, depending on k, but is infinitesimally 
small in magnitude (~ x or v). They are marked with red crosses in Fig. As before, 
the dressed gluon propagator has a branch cut (marked in blue) running from z = —i 
to z = i that is slightly deformed from the imaginary axis. As noted above, in order 
for the analytic continuation to the mass shell limit of the heavy particle to have been 
performed correctly, the z integration contour (green) has to be chosen to go in between 
the two poles of the heavy field propagator and cross the imaginary axis between z = z_ 
and z = i. Outside the vicinity of the origin of the complex plane, we may choose C 
to lie on the real axis. 

For concreteness, let us now define C in the following way: It follows the real axis 
all the way from 2; = — oo to the origin, then — if and only if Z- resides on the upper 
half of the complex plane — makes a detour shooting up along the imaginary axis 
to encircle the pole ai z = z_ in the clockwise direction and returning to the origin, 
and finally proceeds along the real axis to 2; = +00. As the branch cut of the gluon 
propagator has been shifted away from the imaginary axis, the contributions from going 
up and down the imaginary axis cancel, and we are left with two separate parts: One 
from the integral along the real axis, and another from encircling the pole at z = z_ in 
the clockwise direction. We denote these two functions by /^^^^ and /|°^^, respectively, 
and now proceed to expand them in a power series in v while working only to leading 
order in small x. One possible non-triviality in this procedure is related to the order of 
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taking v and x to zero, which we have dealt with by separately considering the cases 
V ~ and v r-' x, and checking that they lead to consistent results. Below we will, 
however, only present the calculations in the former limit, which corresponds to first 
taking the small-x and only then the small-^; limit. 

The pole contribution 

Let us first inspect the pole contribution, Z^"^*^, which is only present for momenta k 
that satisfy v ■ k > xk/2. Denoting the angle between v and A; by ^, we easily obtain 

87r2 7o 7o y/l + (v-xky + f{z-) 

F{x,v), (A.15) 



X 



87r2 

where for v non-zero and of order we may set x to zero in F. This leads to 

F(0, v) = r dk d9 sin g — i /(^^ cos^) 

^ ' Jo Jo VTT^ k^ + fi^v cose)' ^ ' 

in which we expand / in powers of f , using the fact that consistently with our choice 
of the branch cut, we have around the origin 

This enables us to perform the above integrations with ease, giving for the real part of 

jpole 

The real cLxis contribution 

Next, we move on to consider the integral along the real axis. Here, we take advantage 
of the reality of z and write 1^ in the form 



= -X 



dz f d^k 1 f(z) r x(l + z'^)e - 2v k 



_^ 27r J (27r)3 k k"^ + f{z) { (x(l + z^)k^ - 2v ■ kf + 4.z'^k'^{l + v^) 



. 2zkVlT^ 1 

^ ' {X{1 + Z2)A;2 - 2V ■ kf + 4^2^2(1 + ^;2) J' ^ ■ ^> 

where the latter, imaginary part obviously vanishes due to its antisymmetricity in z. 
For the real part, we on the other hand see that for z of order x^, the k integral vanishes 
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by symmetry to leading order in x, so we must let z = xz. Setting then x ^ wherever 
possible, we obtain in terms of this new variable 

r-i = -r r ^ 1 -J^ - 2^ ■ k/x 

^ J_^2nJ{2n)^kk^ + f{0){k^-2vk/xy + 4^k^{l + v^)'^' ' 

This expression is obviously at least 0{x^)^ iiv is of order a;°, and we may hence neglect 
it here. 

Final result 

Assembling the two pieces of our calculation together, we finally obtain the desired 
result 

Re/. ^ Re + H^) ^ "T^^ - X " S ) + °(^''}- f*'^^' 

This result may be verified by considering velocities v of order x, in which case the 
individual results for /^"^'^ and however differ from the above expressions. 

A. 2 The transverse integral 

With the transverse integral. It-, our task is somewhat simpler, as to obtain the rest and 
kinetic masses, it suffices to consider the integral in the zero external three-momentum 
limit, t; = 0. The calculation proceeds in almost exact parallel with the above one, 
with the only subtlety being related to the unphysical Landau pole go = that 
appears when evaluating the transverse propagator on the Euclidean real axis between 
the branch cut and go = 0.^° Due to this pole, one cannot combine the go-integration 
as in going from Fig. 0to Eq. ( |A.6|) , but must evaluate it in two pieces, from go = —oo 
to the branch cut and from go = gL(g) + e to infinity. 

Denoting the integration path described above by C", it is straightforward to see 
that to leading order in x we obtain 



It = 



(P - Q)2 + M2 tg2 + g2 + Ht g^ + g 



IX 



dz 1 [ (Pk I- kl f{z) 



c" 2n z{l + z^) J {27ry' P (1 + ^2)p + j(^) 



(A.22) 



lOc 



Solving for the poles of the transverse gluon propagator between the branch cut and the imaginary 
axis, one finds the so-called Landau pole at go — Qhiq) ~ — ^^^i , where the last approximation 
corresponds to the q ^ limit. 
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where in the vicinity of the origin, z = 0, the function f{z) is given by 

iz(l + z"^) f , 1-iz 



m = ~ - 



4 



log 



l + iz 



(A.23) 



Observing that the integrand in Eq. (|A.22|) is imaginary everywhere on the Euchdean 
real axis, we conclude that all real contributions originate from avoiding the possible 
pole at z = 0. This time it, however, is immediately obvious from the form of f{z) that 
the corresponding residue vanishes (i.e. there in fact is no pole at z = 0), from which 
it follows that 



(A.24) 



A. 3 The self energy 

Assembling everything together, we finally see that we obtain the result of Eq. (|6.43|) 
for the self energy function Efgg(P), evaluated to leading order in x = nici/M and to 
order f ^ in powers of the velocity. 



B. Integral representations of sum-integrals 

Here, we assemble the integral representations of various sum-integrals needed in our 
calculations. Using standard integration tricks, one easily obtains 

1 1 



± (P + K)^ + M2 



dk Re 







log 



°° k 
dk-—-—Re 

E{k) 
1 



log 



P2 + + 2ikpo + 2pk 
P^ + M'^ + 2ikpo - 2pk 
_p2 + ]vf2 + 2iE{k)po + 2pk 
-P2 + M2 + 2iE{k)po - 2pk 



sgn(P) 
emk) ^ sgn(P) 



"27r2 



dk 



.(B.l) 
(B.2) 



where sgn(P) = 1 for bosonic and —1 for fermionic P, E{k) = yP + M^, and the 
operation Re is defined by Re/(po) = 1/2 [f{po) + f{—Po)]- 
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Kinematic region 



M = 0, p < 



M = 0, p > 



gT^M^ T/g, 
p unconstrained 



M > T/g, 

p non-relativistic 

rriD and refer 
to the static 
thermal screening 

masses of the 

gluons and scalars 



Dispersion relation 



QCD: 



2m2p +{£±p)\m'^ logf — - ± 2p'^ 



m 



SUSY: 



quarks: 2m^p + ± p) < log— ± 2p 

c +p 



m 



'S + p 



= 0, 



= 



scalars: = p' + 6M'' 6M^ = g^CpT^ 



QCD: 



SUSY: 



g'CpT' 



QCD: (see Fig. 4) 

£2 =p2^M2 + 5M2, 5M ^5M{p,M) 
SUSY: (see Fig. 2) 

£2 =p2^M2 + 5M2, 5M^5M{M) 



QCD: 



E = Mrest + -Mkin^;^ 



57r 



g^Cpmo 



247r 



TT 

16 



SUSY: 



1 



rest 



g^Cpmp 
247r 



TT 

16 



+ 



Stt 



Table 1: This table summarizes our results. The SUSY results apply for both the funda- 
mental quarks and scalars in the N = 2 hypermultiplet unless otherwise stated. 
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